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Abstract 

We present analytic evidence for the occurrence of an upsilon point, an 

infinite checkerboard structure of modulated phases, in the ground state of 

a spin model. The structure of the upsilon point is studied by calculating 

interface-interface interactions using an expansion in inverse spin anisotropy. 
pacs: 05.50.+q 64.40.Cn 75.10.Hk 
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I. INTRODUCTION 



Simple spin models can have surprisingly complex phase structures, even at zero tem- 
perature. In particular, near multiphase lines, along which the ground state is infinitely 
degenerate, a perturbation such as temperature [ 0], quantum fiuctuations [ or softening 
of the spins [ ^ can result in infinite sequences of stable phases. 

The existence of upsilon-points (T-points) in Frenkel-Kontorova models has been pointed 
out recently [ ^, These occur when two multiphase lines meet at a first-order boundary. 
A small perturbation about such a point can stabilise an infinite checkerboard structure of 
commensurate phases as shown in Fig. |l|. In many ways the T-point can be thought of as 
a two-dimensional generalisation of the behaviour customary near a multiphase point. 

The occurrence of an T-point in a spin model was recently suggested by the numeric 
work of Sasaki [ ^. Here we present the first analytic evidence for the existence of a T-point 
in a spin model. The system we consider is the chiral XY model with 6-fold spin anisotropy 
in a magnetic field [ 0. We identify a candidate for a T-point at infinite spin anisotropy 
D and show that, as D is reduced from infinity, the softening of the potential wells allows 
formation of a T-point structure. 

We follow the method introduced by Fisher and Szpilka [ ^ and extended by Bassler, 
Sasaki, and Griffiths [ |^ and identify the multiphase structures as comprising interfaces 
separating domains of the different phases. The behaviour near a single multiphase point can 
be analysed in terms of a unique type of interface. However there are two (or more) different 
phases stable near an T-point and hence different types of interface must be identified in the 
analysis. It is the interactions between the interfaces which are responsible for breaking the 
multiphase degeneracy and a knowledge of their sign and dependence on separation allows 
determination of the phase diagram. Here the interface-interface interactions are calculated 
using an expansion m D^^ [|^. 
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II. THE MODEL 



We consider the classical chiral XY model with 6-fold spin anisotropy, D, in the presence 
of an external magnetic field h. The Hamiltonian of the system is 

n = J2l- cos{9i^i -9i + 71 A/3) + h[l - cos(^i)] + D[l - cos(6^i)]/36| (2.1) 

i 

where Qi is the angle between the i*^ spin and the magnetic field orientation. We shall 
concentrate on the behaviour of the model near the limit D = oo, where tt,,, defined as 
SOi/iT, can take only the integer values {0,1. ..,5}. 
The ground-state configurations satisfy 

f ^0 «. (2,2) 

For a given i, equation (|2.2| ) enables us to express 9i^i as a function of 9i and 9i_i. This 
fact, together with the observation that rii G {0, l...,5}, is sufficient to conclude that, for 
D = oo, there will always exist periodic minimal energy configurations. It will be convenient 
to label a periodic configuration {..., 9^, 9i, 92, ■■■9n, 9i, ...} as {rii n2 ...un)- 
We can now discuss the phase diagram for D = oo, obtained using the Floria-Griffiths 
algorithm [ |T0|, and presented in Fig. |^. We have restricted the labelling of the phases to 
the first quadrant (0 < A < 3; /i > 0); the remaining phases can be constructed through 
appropriate symmetry operations on the sequences. The transition lines between regions 
A and J and regions J and F are first order. The remaining boundaries are multiphase lines, 
that is loci where all phases (including non-periodic ones) built from arbitrary combinations 
of the two neighbouring phases are degenerate [ [l|]. 

If the spin anisotropy is reduced from oo it seems natural to expect the degeneracy along 
the multiphase lines to be lifted as the spins soften from the clock positions. Although rii 
is no longer constrained to assume integer values, nevertheless, for high values of D, the 
angles 9i will be close enough to the clock positions to allow us to continue to use the same 
labelling scheme. 
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We are particularly interested in the possible appearance of T-points for finite D. An 
T-point can occur when a first-order transition line separating, say, phases (a) and {(3) 
(that, for simplicity, we now assume to be non-degenerate) approaches a commensurate- 
incommensurate transition. An infinite number of phases spring out from the multicritical 
point at the end of the first-order line (as represented in Fig. |I|). The phases appearing are 
made of sequences of (a) and {13). As a 7^ /? the interfaces separating them, which we shall 
call lajs and Ipa are also generally distinct. In the example of Fig. the general form for a 
phase in the fan is {d"'Ia(3f3^Ipa) , where the integers n, m increase approaching the (a) and 
{(3) boundaries respectively. 

The multiphase point P highlighted in Fig. |^ seems to be a good candidate for becoming 
an T-point when D is relaxed from 00. P lies at the end of a first-order transition line and 
it seems reasonable to consider the two multiphase lines J-G and G-F as special cases of 
accumulation lines. Therefore we might expect to observe a structure similar to Fig. |I| for 
small values of 1/D. 

III. THE (a) BOUNDARY 

Consider the J-G boundary (at a finite distance from the point P). When D = 00 the 
phases (a) = (51) and (51402) coexist, and it is easy to check that along the boundary all 
phases built with a sequences separated by a |402| = / block (i.e. (a^/a™...)) are degenerate. 
We want to study how this degeneracy is lifted when D assumes finite values. 

It is physically appealing to regard the / block of spins as an interface separating pure 
a sequences. Indeed one can conveniently write the energy per spin of, say, phase (a"/) as 

E = El + {2n{El - E?) + a + K.(2n) + K.«(2n, 2n) + ...)/(n/ + 2n) (3.1) 

where n/ = 3, E'j {E^ is the energy per spin of phase (/) ((«)), u is the creation energy of /, 
Va{2n) is the interaction energy of two interfaces / separated by a distance 2n, Vaa{2n, 2n) 
is the interaction energy of three interfaces and so forth. In the (a) region the interface 
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tension a is positive; as the phase boundary is approached a decreases and eventually, when 
it is balances the interface interactions, it will be favourable for the system to replace the 
pure a phase with a modulated one. The nature of the transition depends on the form of 
the interface interactions, which we now calculate to leading order using an expansion in 
inverse spin anisotropy. 

In the large D limit Va{2n) dominates the energy contribution from the interface inter- 
action terms. It can be obtained using the reconnection formula [ ^ 

\4(2n) =Ei + E2-E^-E^ (3.2) 

where Ei is the energy of configuration i sketched in Fig. ^. 

Equation (|3.2|) is exact, but is not convenient for our purposes, as we want only the 
leading term of Va{2n). In fact, we can exploit the rapid decay of the Vq, with n to substitute 
all infinite segments in Fig. |^ with finite (though sufficiently long) ones. Thus equation (3.2) 
can be approximated by 

K(2n) ^Ep^ + E^-Ec (3.3) 

where Ai, B, and C, are the periodic configurations sketched in Fig. ^ ni + 2n and — 
2n — 2ni — rii are assumed to be both a large multiple of 2n and much greater than nj. 

We label the spins of configurations A, B and C as {ai, 02, = ao}, {bm+i, = 
bni} and {ci, C2, cat = Cq} respectively. For D large the spins will deviate from their clock 
positions {a^}, {b^} and {c^} by an angle analytic in D^^ and we write 

Oi = a° + di, bi = b°i + bi, q = c° + q. (3.4) 

We can choose to label the spins such that 

= c° 1 < 2 < ni, ^° = c- ni + l<t<N. (3.5) 

Then using a superscript tilde to indicate we are working only to second order in the spin 
deviations {Sj}, {bi}, {cj} the two-interface interaction can be written 
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where 



i=2 «=ni+2 
N _ 

-H{cn, ci) + H{ci-i, Ci) (3.6) 

i=2 



Hiai-i, a,) = Jf_i ,{a,_i - + A^i J' + ^^a,, + e,")^ + D~al/2 (3.7) 



with 



Jti, = cos(a°_i - a° + 7rA/3)/2, (3.8) 

/i^ = /icos(a°)/2, (3.9) 

, = tan(a°_i - a° + vr A/3) , (3.10) 

e^ = tan(a°). (3.11) 



It follows from (|3.5| ) that 



e^ = 4, K = K (3.12) 



for all i and that 



Jti, = ^1,., ^-1, = A^_i„ 2 < ^ < (3.13) 
J' = Jl,„ At, , = Ati„ + 2 < z < iV. (3.14) 



For the cases considered here it will be possible to label the phases so that ( p.l3|) is also true 



for i = 1 and ( 3.14|) for i = ni + 1. Under these circumstances we may drop the a, b and c 



superscripts on the quantities defined in (|3.8| )-( prTl| ). It is then possible to use the recursion 



equations (p.2|) to simplify ( p.6|) . After some algebra one obtains 

Vai2n) = -Jo,i{(ani - ^iv)(ci - c„i+i) - (oi - 6„i+i)(cAr - c„J}. (3.15) 

The quantities appearing in ( ^.15| ) can be obtained to leading order in 1/D using the 
recursion equation ( ^.2] ). An example of how to calculate V^(2n) is given in Appendix A. 
The result for general n is 
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Va{2n) = c^crHs4 - sslV^^" + (3.16) 

where Sj = sin[7r(A — ^)/3] and Cj = cos[7r(A — Terms of higher order than quadratic 

in the Hamiltonian ( |2.1|) will not contribute to the leading term of the interface-interface 
interaction and hence to leading order Va{2n) and Va{2n) will be equal. Therefore we shall 
not distinguish between them below. 

A knowledge of the leading term in the interface-interface interaction, equation ( p. 161) , 



allows us to take the first step in determining the ground-state configurations. Because 
we are considering only two-interface interactions the interfaces must be equispaced in the 
ground state. Inspection of equation ( |3.16| ) shows that Va{2n) is always positive and convex 



near P. This is enough to conclude that, for D large, all transitions (a"/) — ^ occur 
as a is lowered [ §]. 

To this order of approximation the (a""/) : (a""*"^/) phase boundaries remain degener- 
ate and higher-order interface interactions can introduce qualitative changes in the phase 
diagram. This will be discussed further in section V. 



IV. THE (/3) BOUNDARY 

We now focus our attention on what happens along the {(3) boundary in the two-interface 
interaction approximation. In the F region of the phase diagram (Fig. ^) three phases coexist 
when D = oo, namely Pi = (14), P2 = (25) and = (03). However, when D is relaxed, 
only phases /5i and P2 continue to stay degenerate, while phase Ps has a higher energy. 

Consider the boundary between one of the phases (a"/) and region F. Along this bound- 
ary, in the absence of interactions between the interfaces Ji = (51)"4, I2 = and I3 = 2, all 
phases ((51)"4(14)™i0(30)'i2(52)pi(51)M(14)'"20(30)'22(52)P2....) are degenerate. 

Now we turn on the two-interface interactions. In this approximation, the possible ground 
states are periodic and have the form ((51)M(14)™'0(30)'2(52)p), where m, / and p depend 
on cr. In the following analysis we shall hold n fixed and assume that a can be varied to 
trace out the phase sequences. 
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The energy per spin can be written 

^ = I (1 + 2n)Ei^ + Ep^2m + Ei^ + Ep^2p + Ei^ + Ep^2l + a 

+V>,(2m) + Vp,{2l) + Vp,{2p)^/L (4.1) 

where L = (2m + 2p + 2/ + 3 + 2n) and cr includes the energy tension of the three interfaces 
Ji, h and /a. 

Simple calculations show that 

Ep, = Ep„ Ep, = Ep, + 3/iV(8D) + Oil/D-"). (4.2) 

Proceeding as in Section III and Appendix A the two-wall interactions between interfaces 
bounding phases (/32), and {[3^) are to leading order 

Vp,{2m) = Vp,{2p) = Z}-(2"+2){s2 - S3}'cr+\ 

\/ft(2/) = W^,(2/)/c3. (4.3) 



We now want to find the values m, p and / which minimise ( |4.1| ) for a given n and a. 
By symmetry arguments one has m = p. It follows from ( [4.2| ) that / must be bounded 
from above. Indeed an explicit calculation of the energy 0{1/D) shows immediately that 
r = or 1 and that the sequence of phases as a is lowered is, using the notation [n, rh, J, p] , 
[n, 0,0,0] ^ [n, 0,1,0] ^ F. 

The boundary between [n, 0, 0, 0] and [n, 0, 1, 0] is non-degenerate and cannot be split by 
terms of higher order in D~^. The boundary between [n, 0, 1, 0] and F remains degenerate 
to all phases of the form [n, 171,1,771]. The effect of higher order terms can be deduced 
by noting that Vfl^(2m) and Vp^{2p) are positive and convex. This implies that all the 
transitions [n, m, 1, m] [n,m + 1,1, m + 1] are stable [ P]. Fig. 5 summarises the results 
of the two-interface interaction analysis. 



V. HYPERFINE STRUCTURE 



We now restrict our analysis to the richest region of the phase diagram, i.e. where / = 1. 
We already know that, in the two-interface approximation, the possible ground states can 
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be written in the form {a'^Ii(3^l2f3^l3) = [n,m,m], where Ii = 4, J2 = 030, and I3 = 2. 
Bassler, Sasaki and Griffiths [ ^ have shown that for exponentially decaying interactions 
such as is the case here the general form of the interaction energy of an arbitrary number 
of interfaces can be constructed as 

Va/3i/32.../32(2^, 2m, 2p, 2g, . . . , 2s) = 

V^{2n)t^p,Vp,{2m)tp,pM2p)tp,M'^l) " --tp^pM^s) (5.1) 



where the V's are defined in (|3.16|) and ( [4 .31) and to leading order 

ta(3i = D{{S3 - S2)(S4 - S3)}"\ 

^/32« = taPx- (5-2) 

The formulae ( ^.21) follow from calculations similar to that described in Appendix 
A. For example taking the phases A = [n, 171,171], B = [n,m + l,m + 1], and C = 



[n,m,m,n,m + l,m + 1] the right-hand side of equation ( |3.15| ) is equal in leading order 



to Vi3j^{2m)tf3j^/3^V/3^{2m + 2). t^^^j can then immediately be extracted by using the expres- 
sion ( [4. 3D for Va^ and Vf^^. 

With the aid of (|5.1| ) it is possible to examine the effects of three-interface interactions 
on the superdegenerate boundaries in Fig. 5. Consider the general case represented in Fig. 
^(a). All four boundaries are multiphase lines where any sequence of the two neighbouring 
phases are degenerate within the two-interface interaction approximation. For the [n, m, m] : 
[n+1, m, m] and [n, m+1, m+1] : m+1, m+1] boundaries this exhausts the possibilities 

and the three-interface interactions are not of sufficiently long range to split the degeneracy. 
For the [n, m, m] : [n,ni + l,m + l] {[n + l,m,m] : [n + l,m + l,m + 1]) boundary, however, 
the phases [n, m + 1, m] and [n, m,m + 1] {[n + l,m + l,in\ and [n + l,m,m + 1]) are also 
degenerate and there is the possibility that these may be stabilised with respect to [n, m, m] 
and [n, m+1, m+1] ( [n + l,m, m] and [n+1, m + 1, m + 1]) by the three- interface interaction. 

To check this we need the energy differences 
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2-E'[n,m+l,r?i] -^[n,m,m] -^[n,m+l,rra+l] '^'^[n,m,m+l] '^[n,m,m] -^[n,m+l,r?i+l] 

~ Vp^f32{m,m + 1) + Vai/32("^ + l,m) - Val/32("^ + l,m + 1) - Vp^p^{m,m) (5.3) 

which are dominated by Vp^p^ [m, m) and which are therefore negative. Similarly the [n + 
l,m,m] : + l,m + l,m + 1] boundary is unstable with respect to the formation of 
{[n + 1,771 + l,m],[n + l,m,m + 1]}. The resulting modification to the phase diagram is 
sketched in Fig. ^(b). 

The V4 terms do not cause further splitting of the multidegenerate lines of Fig. Kb) but 
they qualitatively change the phase diagram near the two points where four lines meet. In 
proximity of the upper one the structure of the phase diagram is determined by the signs of 
the energy differences [ || 

AVi = Vai3^p^{n,m,m) + Vai3^p^{n + l,m + l,m) - Vafi^p^{n + l,m,m) - Vai3^p^{n,m+ l,m) 
AV'2 = Vp^p^aim, m, n) + Vp^p^a{m + 1, m, n + 1) - Vp^p^aim + 1, m, n) - Vp^p^Jyjn, m, n + 1) 

(5.4) 

With the aid of the factorization formulae ( |5.1| ) it is possible to check that the two energy 
differences ( p.4| ) are positive. This means that phases [^ + 1, vn, m] and {[n, m, m + 1], [n, m + 
1, m]} are separated by a short first-order line; similarly one can show that {[n + 1, m, m + 
1], [n + l,m + l,m]} and [n, m + l,m + 1] also coexist at a first-order transition. In this 
approximation the structure of Fig. ^(b) must be modified as in Fig. P(c). 

The factorization formulae ( ^.1|) allow us to go further and study the the effect on the 
phase diagram of interface-interface interactions of all orders. Bassler, Sasaki and Griffiths 
[ showed that the form of the phase diagram depends upon the sign of the two-interface 
interactions ( p.l6| ) and ( |4.3| ) and the t's, equation ( |5.2| ). Here these are all positive corre- 



sponding to a case where the superdegenerate boundaries at the end of the first-order lines 
in Fig. 6(c) split under the effect of higher-order interface-interface interactions, giving rise 
to a structure analogous to that in Fig. 5 (but where the phases have longer periodicity). 
Furthermore one can carry the analysis further by studying again the splitting near the 

10 



points where four lines meet and so on, finding a structure similar to the one in Fig. 6(c). 
The analysis can then be repeated ad infinitum, showing that the T-point has, indeed, a 
self-similar, fractal structure. 

VI. DISCUSSION 

The analysis presented above was based on retaining only the leading order term in the 
interface-interface interaction. We cannot rule out the possibility that the neglected higher- 
order contributions could affect the phase diagram. In particular there will be correction 
terms 0{P/D^) where / is the period of a given phase which could introduce qualitative 
changes for I large and D not sufficiently small. 

The results were checked numerically in two ways. Firstly we used the Floria-Griffiths 
algorithm [ |T^ on a grid of size 1200 to check which phases appeared. It was possible 
to resolve phases with I up to 13. Secondly we used a mean-field analysis, exact at zero 
temperature, to check the positions of the phase boundaries. In this way the the formulae 
for the interface-interface interactions could be verified for short-period phases (typically I 
up to 9). 

To summarise, we have presented analytic evidence for the existence of an T-point in 
a spin model. The phase diagram has been constructed inductively by calculating the 
interface-interface interactions to leading order in 1/D, the inverse spin anisotropy. Follow- 
ing arguments due to Bassler, Sasaki and Griffiths [ |] we have argued that the T-point has 
a self-similar, fractal structure. 
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APPENDIX A: CALCULATION OF THE TWO-INTERFACE INTERACTION. 



As an example of how to obtain the two-interface interaction we consider exphcitly the 
calculation of Va{Q)- Following Fig. 4 we need to consider the periodic phases listed below 
where rii = 4, N = 24 and n = nj = 3. A choice of labelling that satisfies ( ^.131 ) for 
1 < i < rii and ( p.l4|) for rii + 1 < i < is as shown. 



A: 3a?/7r 



i = ni; 1, 2 ... 

i i 

5 15 1 



B: 3b°/7r 



... N;ni + 1 ... 
i i 

5 15 15 140251515151402 



C: 3c0/7r 



N; 1 ... ni,ni + 1 ... 

J, J, J' J' 

5 15 15 15 15 140251515151402 



(Al) 



We can now use ( p.l5| ) to calculate Va{6) to leading order. The quantities (a„^ — b^), 
(5,1 — bni+i), (ci — c„j+i) and (cat — c„J can be obtained correct to leading order from the 
linear approximation to the recursion equation 



where we have used the definitions (|3.8|)-( prTTD 
Let 



(A3) 



Substituting into (|A2|) and equating like powers of D ^ gives 
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To calculate (Si — &ni+i) it is helpful to display explicitly a° and &° as a function of i. 



i ... -3 -2 -1 1 2 3 4 ... 

3a?/7r ... 5 1 5 1 5:15151 5 1 ... 
3^ni+i/^ •••0 2 5 1 5:1 5 1 4 2 5 ... 

(A6) 

The *'s mark where a° first differs from when moving away from the dotted interface 
in either direction. 

It follows immediately from (^5) that 



(S3 - K+s) = ^{mI^ - 3] - sin[|A - 4]} (AT) 



Two further iterations of (^5) give 



~ 1 TT TT ~ 

(ai - = — cos[-A + a? - a"] cos[-A + a" - al]{d3 - 



= cos[^(A - 4)] cos[|(A - 2)]{sin[^A - 3] - sin[^A - 4]}. (A8) 
(a„j — &Ar) may be calculated in an analogous way 

(Sni - &Af) = -(Si - ^ni+l)^ (A9) 

Similarly 

(ci - c„^+i) = (cat - c„J = (ai - (AlO) 

Using Jo,i = cos[|(A — 2)]/2, from the definition (|3l8| ) and substituting (|A8|) -( |AT0|) into 
( p.l5| ) gives 
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\4(6) = -^{sm[^A - 4] - sin[^A - 3]}' cos[^(A - 4)]^ cos[^(A - 2)]\ (All) 

It is important to point out that the labelhng used in equation ( [Al| ) is not unique. Any 
labelhng which satisfies ( p.l3| ) for 1 < i < ni and ( |3.14| ) for 1 < rii + 1 < will give the 



correct results for V^(6). However in general (ai — a„^+i) etc. will contain lower powers 
oi 1/D which cancel when the difference in ( p.l5| ) is taken. The choice given above, which 



maximises the distance of the position (*) where a^ first differs from b^^_|.^ avoids such 
cancellation and leads to the easiest calculation. 

Finally we mention that because the interface-interface interactions decay very rapidly 
(exponentially) with increasing interface-interface distance the values of ni + 2n and 
N — 2n — 2ni — ni need not in fact be much larger that 2n and rij. That sufficiently 
large values have been chosen can be checked a posteriori by verifying that increasing the 
values of ni and A^ does not change the result ( |A11| ). 
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FIGURES 
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FIG. 1. Schematic representation of an T-point. 




FIG. 2. Ground state of the Hamiltonian (2.1) for D = oo. A = (0); B = (012345); 
C = {(024) ,(153)}; D = (01245); E = (0135); F = {(03) ,(14) ,(25)}; G = (02514); H 
= {(025), (014)}; I = (015); J = (15) 



1) 
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3) 



4) 



FIG. 3. Configurations needed to calculate the two-interface interaction; see equation (3.2). 
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FIG. 4. Periodic configurations appearing in the approximate reconnection formula (3.3). 




FIG. 5. Schematic phase diagram near the T-point in the two-interface interaction approxima- 
tion. The notation [n,m,l,p\ is used to indicate the phase ((51)"4(14)"*0(30)'2(52)p). The bold 
solid lines are accumulation lines. The dashed line is a first order boundary. 
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FIG. 6. Detail of the phase diagram in (a) the two-interface approximation; (b) the 
three-interface approximation; (c) the four-interface approximation. First order hnes are dashed. 
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